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Abstract 

Let u be a non-negative super-solution to a 1-dimensional singular 
parabolic equation of p-Laplacian type (1 < p < 2). If u is bounded 
below on a time-segment {y} x (0, T] by a positive number M, then it 
has a power-like decay of order with respect to the space variable x 
in R X [T/2,T]. This fact, stated quantitatively in Proposition 11.11 is a 
“sidewise spreading of positivity” of solutions to such singular equations, 
and can be considered as a form of Harnack inequality. The proof of such 
an effect is based on geometrical ideas. 
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1 Introduction 


Let E = {a,/3) and define E-t„,t = E x {—To,T], for Tq, T > 0. Consider the 
non-linear diffusion equation 

= 0, l<p<2. (1.1) 


A function 


u G Cioc( - To, T; LL(L;)) n LL( - r„, T; (A)) (1.2) 

is a local, weak super-solution to (HU, if for every compact set K C E and 
every sub-interval [ti,t 2 ] C (—To,T] 





Px] dxdt > 0 


(1.3) 


for all non-negative test functions 

e ( - To,T-L\K)) n LL( - r„,T; Wl^^{K)). 

This guarantees that all the integrals in (11.31) are convergent. These equations 
are termed singular since, for 1 < p < 2, the modulus of ellipticity \ux\^~^ —t oo 
as \ux\ 0. 


Remark 1.1 Since we are working with local solutions, we consider the domain 
E-r„,T = Ex {—To, T] , instead of dealing with the more natural Et = A x (0, T], 
in order to avoid problems with the initial conditions. The only role played by 
To > 0 is precisely to get rid of any difficulty at t = 0, and its precise value plays 
no role in the argument to follow. 

Proposition 1.1 Let u be a non-negative, local, weak super-solution to (HU 
in E-xo,t, in the sense of (fm-IOl). satisfying 

u{y,t)>M VtG(0, (1.4) 

for some y € E, and for some M > 0. Let p = I ) , take p > 4p, and 

assume that 


Bp{x) C B 4 p{y) C E, where dist(a;,y) = 2p. 

There exists a G (0,1), that can he determined a priori, quantitatively only in 
terms of the data, and independent of M and T, such that 


u{x, t) > aM 



T T 

for all {x,t) G Bp{x) x 


(1.5) 
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1.1 Novelty and Significance 

The measure theoretical information on the “positivity set” in {y} x (0, 
implies that such a positivity set actually “expands” sidewise in R x 
with a power-like decay of order with respect to the space variable x. 
Although considered a sort of natural fact, to our knowledge this result has 
never been proven before; it is the analogue of the power-like decay of order 
with respect to the time variable t, known in the degenerate setting p > 2 

(see [2], [SI Chapter 4, Section 4], |7]). As the -decay is at the heart of 
the Harnack estimate for p > 2, so Proposition 11.11 could be used to give a more 
streamlined proof of the Harnack inequality in the singular, super-critical range 
< P <2. This will be the object of future work, where we plan to address 
the general A^-dimensional case. 

The proof is based on geometrical ideas, originally introduced in two dif¬ 
ferent contexts: the energy estimates of § [2] and the decay of § [3] rely on a 
method introduced in |8] in order to prove the Holder continuity of solutions to 
an anisotropic elliptic equation, and further developed in ISIIS]; the change of 
variable used in the actual proof of Proposition 11.11 was used in |3]. 

1.2 Further Generalization 

Consider partial differential equations of the form 


ut - {A{x, t, u, Ux))x = 0 weakly in 


( 1 . 6 ) 


where the function A : x M x R —>■ M is only assumed to be measurable 

and subject to the structure condition 



(1.7) 



where 1 < p < 2, Co and Ci are given positive constants. It is not hard to show 
that Proposition ll.ll holds also for weak super-solutions to (HHi-dni) , since our 
proof is entirely based on the structural properties of (di), and the explicit 
dependence on Ux plays no role. However, to keep the exposition simple, we 
have limited ourselves to the prototype case. 

2 Energy Estimates 

Let u be a non-negative bounded, weak super-solution in assume 


0 < pL- = inf u, 


and let a; be a positive parameter. Without loss of generality we may assume 
that 0 £ (q;,/ 3). For p sufficiently small, so that (—p, p) C {a,l3), let 


Bp = {-p,p), Q = BpX {0,T], 
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Bp{y) = {y- P,y + p), Q{y) = Bp{y) x (0,T], 
a G (0,1), H G (0,1] parameters that will be fixed in the following, 
A = {(x,t) G Q{y) : u{x, t) < p- + {1 — a)Hu}}. 


Proposition 2.1 Letu be a non-negative, local, weak super-solution to CH) in 
E-ro,T, iTT- the sense of There exists a positive constant 7 = 'yip), 

such that for every cylinder Q{y) = Bpiy) x (0, T] C and every piecewise 

smooth, cutoff function ^ vanishing on dBpiy), such that 0 < C < Ij and 0, 


I Bp{y)n{Ax{0}) 


{u — p- -\- ojjjH)'^ p 


2-p 


iujH)P-^ 


C^ix, 0 ) dx 


Q{y)nA 


(it — /i_ + aujH)P 


(P dxdt < j [[ \C,xy’dxdt ( 2 . 1 ) 
JjQ{y}nA 


7 / / {u — p--\-aujH)'^ P^P ^\Q\dxdt. 

'Qiv)nA 


Proof - Without loss of generality, we may assume y = 0. In the weak formula¬ 
tion of (ED take p = G{u)(P as test function, with 


G(u) = 


(u — p- -\- aujH)P~^ iuiH)P~'^ 


and C = Ciix)f, 2 it)^ where Ci vanishes outside Bp and satisfies 

0<Ci<l, Ci = linB£, |9:.Ci|<—, 

" P 


( 2 . 2 ) 


for an absolute constant 71 independent of p, and ^2 is monotone decreasing, 
and satisfies 

T 'lo 

0 < C 2 < 1, C 2 = 1 in ( 0 , -], C 2 = 0 for t > T, | 5 tC 2 | < f , (2.3) 

for an absolute constant 72 independent of T. It is easy to see that we have 

p-\ 


= -7 ^ mP 

[u — p- -\- aujH)P 

Modulo a Steklov averaging process, we have 


XA- 



UtG{u)(P dxdt 


JJ (PG'{u)\ux\^ dxdt-{-p JJ Giu)\ 


|P-2 AP-l 


C Ux • Cxdxdt > 0 , 
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(p-1) 


QnA (w - P- + au}H)P 


dxdt 


<pII C^-'t-^ 

iQnA [u-fi-+au} 


Ut 


QnA (w - P- + au}Hy 


j^lUdxdt 

—C dxdt - [[ mp-i C^ dxdt, 

^ JjQnA i^H)P 1 


(P- 1) 


I QnA (m - p- + au}H)P 

<p 11 

IQnA 


dxdt 


dt 


QnA 


/ I ttNd-I I 

(u — + aujH)P ^ 

(u — ^- + au}H)'^~P u 


(P- 1) 

< P 


2-p 


QnA {u- ^l-+ au}H)P 

^p-1 

QnA 


(wiP)P-i 

dxdt 


dxdt, 


+ 


L 


Bf,n{Ax{T}) 

f 

Bpn(Ax{0}) 


f , dxdt 

[u — fi- + aujH)P ^ 

{u — ^- + aojHY~P u 


2—p {ujH)p~^ 

(m — /i_ + au;H)'^~P u 


2-P 


CP{x,T)dx 
CP{x, 0) dx 


■P 


[f 

{u — P- + QUjHY P 

u 

JjQnA 

1 

to 

1 

{ujH)p-\ 


{wH)P-^ 

dxdt. 


The second term on the right-hand side vanishes, as ({x,T) = 0. Therefore, an 
application of Young’s inequality yields 


(p-1) 


QnA {u- H-+ aujH)P 

{u — p- + aujH)‘^~P 


(P dxdt 


< 


JBpn(Ax{0}) 
p-1 


2-p 


(P^x, 0) dx 


QnA (m - P- + aujH} 


{ujH)p-\ 

-CPdxdt + ^ (( \C,x\P dxdt 
JjQnA 


I QnA 2-p 
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where we have taken into account that Q < 0. Therefore, we conclude 


+ 


Bpn(Ax{0}) 

p-1 


[u — fi- -\- aojH)^ ^ 


2-p 


0 ) dx 


2 J JQnAiu- p-+ awH)‘ 


-<^^dxdt<j [( \C,x\^ dxdt (2.4) 

JJonA 


+ 7// {u — p- + oloH)'^ ^\Ct\dxdt, 

J JQnA 

where the first term on the left-hand side is non-negative, since for 0 < s < 
1 — a < 1 the function /(s) = - is monotone decreasing, and /(I — a) = 

> 1 . ■ 


Remark 2.1 The constant 7 deteriorates, as p 1. 


Remark 2.2 Even though in the next Section H basically plays no role, we 
chose to state the previous Proposition with an explicit dependence also on H 
for future applications. 


3 A Decay Lemma 

Without loss of generality, we may assume p- = 0. Let M = w, L < and 
suppose that 

u{0,t)>M Vte(0,^]. (3.1) 

Now, let So be an integer to be chosen, and define 

= {t€ (0,j] : 3 xGBe, u{x, t) <^} 

A{t) = {x€Be : u{x, t) < L(1 - ^)}, t e (0, 
A=Ux,t)&B^ X (0,^] : u{x,t) <T(1-^)|. 


Lemma 3.1 Let u be a non-negative, local, weak super-solution to dnD in 
E-xc,t, in the sense of (11.21) - (11.31) . Let (13.11) hold and take 


L < min{ —, — }. 


Then, for any v G (0,1), there exists a positive integer So such that 
\{tG (0,^] : 3x€Be, u{x,t) < ^}| <i/|(0, ^]|. 
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Proof - Take t £ Ag^-. by definition, there exists x £ such that u{x,t) < 
L/2®°. On the other hand, by assumption u{0,t) > 2T, and therefore, u(0,t) + 
(L/2"->) > L. Hence 


u(0,t) + oST 

ln+ —-^> (so - 1) In 2, 


i{x,t) + 2 ^ 


and we obtain 

(so — 1) In 2 < ln+ 


u{xA) + 


2^0 


d 


= ^ In 


— Inj 


L 


< 


dx 


d 


In 


( L 


^U( 0 ,t) + 

]) df 


dx \ \ u{x, t) + 


dx 


lBnnA{t) 




L 


l{x,t) + 2 ^ 


dx. 


If we integrate with respect to time over the set Ag ^, we have 


(So 


l)|H,Jln2 < 


/o JB^nAit) 




l{x,t) + 2 ^ 


dx 


< 


'0 JBnnA{t) 




L 


l{x, t) + ^ 


dxdt 


< 


QnA (m + 


-fP dxdt 


iQr 


_i 

p 



Apply estimates (12.11) with a = Huj = HM = L. The requirement H < 1 
is satisfied, since L < ^. They yield 


(So 


l)|AsJ < qlQr*’ 
+ l\Q\^ 


QnA 


\C,xY‘dxdt 


ff ^\Ct\dxdt 

J JQnA ^ 


1 

p 


With these choices, we have 


(So 


l)|AsJ < qlQr*” -\Q\^^ +l\Q\’’^ 

P 


L2-P 

T 


_1 

P 


IQI 


_i 

p 


7 




























< 7 


1 fL'^-P 
- + I -jr- 


IQI- 


If we require L < 
we have 


T \ 


, and we substitute it back in the previous estimate, 


T, 


(so-l)|A«„ <7i|(0, ^]|- 

T 

Therefore, if we want that |^s | < ^^1(0, —]|, it is enough to require that So = 

^ + 1 . 

V 

The previous result can also be rewritten as 

Lemma 3.2 Let u be a non-negative, local, weak super-solution to (ED in 
E-ro,T, i‘^ the sense of (jl.2p - f|1.3p . Let (EH hold. For any v G (0,1), there 
exists a positive integer So such that 


T 


|{t e (0,-] : 3a; e Be, u{x,t) < — < J^|(0,- 


J' \ 2^p 2 


provided p > 0 is so large that — ^ TT- 


Now let p be such that 


T\2-p M 


= 


/'2^~pT\ p 


(3.2) 


\M^- 

and assume that Bp C {a, (3). Then Lemmas l3.1113.21 can be rephrased 

Lemma 3.3 Let u be a non-negative, local, weak super-solution to (ED in 
E-ra,T, in the sense of (11.21) - (II.3p . Let (13.11) hold. For any v G (0,1), there 
exists a positive integer So such that for any p > p 


T M f p\~ T 

|{tG(0,-]: 3xGBe, u{x,t)<^[f^ }|<^|(0,-]|, 

provided that Bp C {a,/3). 

Remark 3.1 The previous corollary gives us the power-like decay, required in 
Proposition 11.11 

Let us now set (0, Then, if (13.11) holds, we conclude that for 

any t G and for any x G Be with p > p 


u{x,t) > 


M f p\ 


2 So+i yp 


8 


(3.3) 








Let c > 4 denote a positive parameter, choose x G {a, 13) such that |a;| = 2cp, 
and consider Bcp{x). Then, by (lO) 

M / 2 \ ^ 

provided (EU holds, and B[x)cp C (a,/?)- 

4 A DeGiorgi-Type Lemma 

Assume that some information is available on the “initial data” relative to the 
cylinder B 2 p{y) x (s, s + 9f]P], say for example 

u{x,s)>M for a.e. x € B 2 p{y) (4.1) 

for some M > 0. Then, the following Proposition is proved in [5J Chapter 3, 
Lemma 4.1]. 

Lemma 4.1 Let u be a non-negative, local, weak super-solution to uni), and 
M be a positive number such that (HU) holds. Then 

u> a.e. in Bp{y) x (s, s + 0(2p)P], 

where 

e = 5M^-P (4.2) 

for a constant S G (0,1) depending only upon p, and independent of M and p. 

Remark 4.1 Lemma l4Tl is based on the energy estimates and Proposition 3.1 
of [T], Chapter I which continue to hold in a stable manner for p 1. These 
results are therefore valid for all p > 1, including a seamless transition from the 
singular range p < 2 to the degenerate range p > 2. 

5 Proof of Proposition 11.11 

Fix y € E, define p as in (13.21) . and choose a positive parameter C > 4, such 
that the cylindrical domain 

B2^ {y)x{t),^]cE.r„,T. (5.1) 

2 p Cp Z 

This is an assumption both on the size of the reference ball B 2 ^ (y) and 

2 cp 

on T: nevertheless, we can always assume it without loss of generality. Indeed, 
if dSH) were not satished, we would decompose the interval (0, in smaller 
subintervals, each of width r, such that dED is satished working with p replaced 

by 

2^-Pt 
M^-P 
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The only role of C is in determining a sufficiently large reference domain 

5 iv) C E, 

2 p Cp 

which contains the smaller ball we will actually work with, and will play no 
other role; in particular the structural constants will not depend on C. 

Now, introduce the change of variables and the new unknown function 


z = 2 p -, 


_e-^ = -_ 2 . 

C/ rp . 

T 


v{z,t) = ^u{x,t)e^-p . (5.2) 


This maps the cylinder in (15.11) into Be x (0, oo) and transforms (11.11) into 


Ur - (K 


IP-2, 


5)2 — 


1 


2-p 


weakly in Be x (0,oo). 


The assumption dn of Proposition 11.11 becomes 

' 1 '( 0 ,t) > for all r G (0,+ 00 ). 

Let To > 0 to be chosen and set 

k = . 


(5.3) 

(5.4) 


With this symbolism, (15.41) implies 

i'(0,t) > k for all r G (to,+oo). (5-5) 

Now consider the segment 

I ‘^^{0} X {to,To + 

Let u' = \ and So be the corresponding quantity introduced in Lemma l3.II We 
can then apply Lemmas 13.1H3.31 with T = k‘^~P, M substituted by fc, 

1 k 

As, = {tG{to,To + -]:3zGBp,v{z,t)<^^} for p > p, 
ci0f 2—p 

with p = 2 p . Therefore, if c > 4 denotes a positive parameter, we choose 
z G Be such that \z\ = 2cp, and consider Bcp{z), by (j3.3|) 

yzGB^liz), ^tgAI v{z,t)>^^(^^ , (5.6) 

provided B{z)cp C Be- Summarising, there exists at least a time level ri in the 
range 

To <Ti <To + (5-7) 

such that 

'^zGB^i{z), v{z,Ti)>a 


nC'^-P 


where 


1 


2 ®°+i 


Remark 5.1 Notice that CTo is determined only in terms of the data and is 
independent of the parameter Tq, which is still to be chosen. 
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5.1 Returning to the Original Coordinates 


In terms of the original coordinates and the original function u{x, t), this implies 

■’"I ~'^o Hpf 

> <7oMe~ 2-p = Mo in B^p{x) 

where the time ti corresponding to ri is computed from (15.21) and dSH), and 
dist(x,?/) = 2cp. Now, apply Lemma l4.ll with M replaced by Mo over the 
cylinder B^^{x) x + 9{cpY\. By choosing 

0 = SM^~P where S = (5(data), 

the assumption (14.21) is satisfied, and Lemma 14 .1 1 yields 


> hMo = \aoMe 


> 


1 


— o—ze 


for all times 


2«->+2 V 5 c 
ti<t<ti+d 


e 2-p 


M 


in Bop{x) 


1 


2®°(2-p) V5 

Notice that (15.81) can be rewritten as 


-{Pl-Pa): 


with 


> cr ( - 
\P 


_ def 1 
a = 


M in Bp{x), 


2®°+2 V 5 


2\ 


e 2-p' 


(5.8) 

(5.9) 

(5.10) 

(5.11) 

(5.12) 


If the right hand side of (15.91) equals then (15.81) holds for all times in 
— £M‘^~p{cpY , where e = ‘ 

taking into account the expression for p and cto, we conclude that (15.81) holds 
for all times in 

<-) 

Thus, the conclusion of Proposition fLll holds. provided the upper time level 
in (15.91) equals ^• The transformed Tq level is still undetermined, and it will be 
so chosen as to verify such a requirement. Precisely, taking into account (lOl) 

2«°(2-p) 15/ 2 V2, 


T T 

= -(ti - -) = 5 


This determines quantitatively Tq = ro(data), and inserting such a To on the 
right-hand side of (15.111) and (15.131) . yields a bound below that depends only on 
the data; (15.111) and (15.131) have been obtained relying on the bound below for 
u along the segment {y} x (0, T], However, the same argument on the bound 
along the shorter segment {y} x (0, s] for any ^ < s < T yields the same result 
with T substituted by s: the proof of Proposition II.II is then completed. ■ 
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5.2 A Remark about the Limit as p —)■ 2 

The change of variables (ICT) and the subsequent arguments, yield constants 
that deteriorate as p —2. This is no surprise, as the decay of solutions to 
linear parabolic equations is not power-like, but rather exponential-like, as in 
the fundamental solution of the heat equation. 

Nevertheless, our estimates can be stabilised, in order to recover the correct 
exponential decay in the p = 2 limit. However, this would require a careful 
tracing of all the functional dependencies in our estimates, and we postpone it 
to a future work. 
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